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ABSTRACT 

In this paper uniformly rotating relativistic rings are investigated analytically utilizing 
two different approximations simultaneously: (1) an expansion about the thin ring 
limit (the cross-section is small compared with the size of the whole ring) (2) post- 
Newtonian expansions. The analytic results for rings are compared with numerical 
solutions. 
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INTRODUCTION 



. The problem of self-gravitating, axially symmetric rings in 
' equilibrium can be tackled in Newtonian gravity by expand- 
, ing it about the thin ring limit, where the cross-section of the 
I ring tends to a circle. Doing so, Kowalewsky (1885), Poincare 

(1885a,b,c) and Dyson (1892, 1893) obtained series for ho- 
, mogeneous rings and Ostriker (1964a,b, 1965) and Petroff 
' & Horatschek (2008) for polytropic rings. Alternatively by 

using a Roche model (Roche 1873) rings with a sufficiently 

soft equation of state can be described approximately, see 

Horatschek & Petroff (2009). 

With the help of numerical methods it was possible 

■ to study homogeneous rings and their connection to the 

■ Maclaurin spheroids (Wong 1974; Eriguchi & Sugimoto 
, 1981; Eriguchi & Hachisu 1985; Ansorg, Kleinwachter & 
' Meinel 2003c) as well as non-homogeneous rings (Hachisu 
' 1986). Furthermore relativistic rings and their transition 

to the extreme Kerr Black Hole (Ansorg, Kleinwachter 
& Meinel 2003b; Fischer, Horatschek & Ansorg 2005; 
Labranche, Petroff & Ansorg 2007) can be calculated to 
near-machine accuracy. 

However no analytic work has been done for relativistic 
rings. The reason may be that because of the non-linearity of 
Einstein's field equations the expansion about the thin ring 
limit does not work. In contrast, if one additionally expands 
the rings in a post-Newtonian series, following the meth- 
ods used for studying post-Newtonian Maclaurin spheroids 
(Chandrasekhar 1967; Bardeen 1971; Petroff 2003), it does. 
In this paper the first post-Newtonian corrections for homo- 
geneous rings are calculated and the results are compared 
to numerical ones generated using a version of the code 
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described in Ansorg et al. (2003a), but modified to rings, 
cf. Meinel et al. (2008). 



2 BASIC EQUATIONS 

The matter model is that of a perfect fiuid, i.e. the energy- 
momentum tensor reads 



rife ik , I , V \ i h 

= pg + {fi+ — ]u u , 



(1) 



where p is the pressure, c the speed of light, /ic^ the en- 
ergy density and u' the four-velocity. We are considering 
homogeneous matter, by which we mean the mass-density is 
constant; 



jj, = constant. 



(2) 



An axially symmetric, stationary spacetime containing 
(only) a rigidly rotating fiuid can be described using the line 
element 



ds' 



e"(d£i + dz ) + Q e '' [dip — cu dt) 



(3) 



where the metric functions depend only on g and z. The 
coordinates are uniquely fixed by requiring that the metric 
functions and their first derivatives be continuous at the 
fluid surface. We introduce the further coordinates r and x 
(see Fig. 1) 



Q = b — r cos X, z — r sin x, 



(4) 



where b is the 'centre of mass' of the ring's cross-section, 
defined by 



b-- 



Jj ngdgdz 
JJ fidgdz 



(5) 



and we have retained fj, in this and the next equation in this 
paper to have expressions valid for the more general case in 
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Figure 1. A sketch providing the meaning of the coordinates 



which /I is not constant. The parameter 6 is a coordinate de- 
pendent quantity and has no particular physical significance 
in contrast to invariants such as mass or angular momentum. 
Inserting (4) in the numerator of (5), one finds that 



fir cos X dr dx ~ 



(6) 



holds. Since we are considering rigid rotation, we can in- 
troduce the constant angular velocity fl = u'^/u*. A locally 
non-rotating observer would measure the following three- 
velocity for a fluid element 



V := g{fl — ijj)e'^ 



(7) 



see Bardeen (1970). 

It is convenient to introduce the following dimensionless 
quantities 



V- -, u ■- —, il ■- — , 
c il c 



P--=^- (8) 



jic^ 



The field equations are 
47rG/ie2" 



Av 



1 + v'^ 



1 — fi^ 



(l+p)-K2p 



Vu-V{k + u) + ^^e ^ ' (Vw) , 



Ax{k + v) = 



167rG^ . 2a 

5 — PC 



[V(k + , 



A2CJ — 



A_ia 



167rGM(l~c:')(l+p) ^2. 
- Vii- V(3k-!/), 






(9) 



(10) 



(11) 



(12) 



where G is the gravitational constant. Here V is the nabla 
operator in a three-space with the cylindrical coordinates 
{g, if, z) and A^ is a generalized Laplacian operator, which 
when applied to an axially symmetric function / = f^g, z) 
gives 

9V ''"- + 1 df d^f 



^-f-9^ + 



+ 7j4- (13) 



dg dz^ 



The operator is linear and becomes the standard Laplacian 
for m = 0, i.e. Ao — A. The vanishing divergence of the 
energy-momentum tensor, T' ;fe — 0, gives 

{l+p)c''^l-v^ = 1-7, (14) 



where 7 is a constant related to the relative redshift Zq via 



Zo = 



7 



1-7 



(15) 



To be more precise, this is the redshift of photons emitted 
from the ring's surface and observed at infinity that carry 
no angular momentum. 



3 POST-NEWTONIAN EXPANSION 

Let us introduce the relativistic parameter 
._ VGjlb 





C 


l^-Lu; 


and the following expansions 








00 

1=1 


(17) 


00 
\ ^ 2i 

a = 2_, <^2i£ , 


00 
V^ - 2i 

U — y ^ U)2i£ , 


(18) 


00 


00 
1=1 


(19) 


00 
\ ^ - 2i-l 
V = } ^V2i-l£ , 


00 
\ ^ - 2i 

P=2^P2i£ , 


(20) 



where the even/odd powers of e that appear in the sums 
come about from stationarity (simultaneous reversal in the 
direction of time di — >■ — di and in the direction of rotation 
Q — >■ —Q). The surface of the fluid, deflned by vanishing 
pressure, will be parameterized by the function r — rs(x), 
which we also expand 



s(x)=E^2,(x)e' 



(21) 



When one replaces Gji/c^ in the field equations by e'^/fo^ 
and expands them in e, then the leading order equation 
comes from (9) and is the familiar Laplace/Poisson equa- 
tion (in a somewhat unfamiliar notation) 



Au^ = 0, A4 



An 
fe2"' 



(22) 



and the expansion of (14) yields the Bernoulli equation 



P2 



vi 



^2 = -72, 



(23) 



where v\ — gCli/b results from (7) and the superscripts 'i' 
and 'o' refer to the regions inside and outside the fiuid re- 
spectively. To leading order, (6) becomes 



n)(x)cosxdx = 0. 



(24) 



On the Newtonian boundary ro(x), P2 = holds and the 
function !/2 and its first derivatives must be continuous. It 
follows from asymptotic fiatness that 1^2 tends to zero at 
infinity. No analytic solution for rings is known oven in the 
Newtonian case (described by the above equations) . Instead, 
we solve them using an expansion about the thin ring limit as 
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in Petroff & Horatschek (2008). For example, the expansion 
for the Newtonian surface function is 



n)( 



(x) = a 1 + 5Z 5Z Afc cos(A:x)cr* + o{a 
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(41) 
= 0. (42) 



i/i - 1^1 ) =0, 

r = ro(x) 



d{i^i - 1^1) 



with 



and where it will be convenient to introduce 
A 



ln^-2 



(25) 



(26) 



(27) 



dr 



r2 



d'{4^u^) 



dr^ 



for later. The expression for ro makes use of the fact that 
the cross-section tends to a circle of radius a for a — > 0, 
as mentioned in the introduction. The solution given for 
Newtonian rings in the above paper can be used here, where 
the relations 



2 2 

c e V2, 



u 



Pn = fj,c^e^p2, 



On = cb eQi, 
iiN = ceiii 



Vq 



2 2 

-c e 72, 



(28) 



have to be taken into account and the index 'N' denotes a 
Newtonian quantity here to distinguish it from the relativis- 
tic quantity. 

Turning now to the first post-Newtonian correction, we 
find the equations 



Ai(k2+Z^2) =0, 

A-i{a2 + u2) = 0, 

A2(i2 = 0, 

A2CJ2 



167r 



A:y2 = 0, 



A 4 



8n 

62 



-1^2 + ^1 



TP2 



The first two of these equations are solved by 

K2 = ^2 = —1^2, 



(29) 
(30) 
(31) 

(32) 

(33) 

(34) 

(35) 



where we remind the reader that 1/2, P2 and vi are already 
known from the Newtonian order (cf. (28)) and in (34) 02 
was replaced by — !^2- Equation (14) gives 



P4+P2 \ V2 



vl 



■Vi 



+ 



V2 



1 i-2 Vt 

2-2«i - ^ 



(36) 



-74, 



where 



W3 



f23-f2i((I;2 + 24) 



(37) 



follows from the expansion of (7). Equation (6) yields 



r2{x)n){x) cosxdx = 0. 



(38) 



Additionally, the transition conditions at the fluid boundary 
become 



aj2 - cDa ) =0, 



d{Cb2 - UI2) 

dr 



0, 



(39) 
(40) 



r = ro{x) 
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r = ro(x) 



r^raix) 

at this order. Vanishing pressure on the surface translates 
to 



P4|r=ro(x) +r2{x) 



dp2 
dr 



0. 



(43) 



r=ro{x) 



4 EXPANSION ABOUT THE THIN RING 
LIMIT 

As with the Newtonian order, we expand the unknown func- 
tions as a power series in a and a Fourier series in x (the 
sine terms do not appear as a result of reflectional symmetry 
with respect to the equatorial plane): 



^4 = ^'^' E E ^^^(y) cos(fcx)'^' + o{a'+^), (44) 

^''^'E E F,My,x)a''-'a^-'+o{a''+''), (45) 

(46) 
^^2 = ^E E GuKi{y,x)a''+'a^-'+o{a''+^) (47) 

1 = 1 i=l-l 
q i 

PA = 7r^o-''EE*'=(2/)cos(fcx)o-'-|-o(a'+''), (48) 

i=0 fe=0 

~ n''/^a^Y.Ua^ + o{a''+''), (49) 

1=1 

74 = 7r2a"E5^^' + «(^'^')' (50) 



V4, 



CJ2 



i=0 fe=0 

9+1 q 

1 = 1 i = l-l 
q i 

TTo-^ E E ^'^kiv) cos{kxW + o{o 

<?+i q 

21+1 

-^u^i(y,XK 

1 = 1 i=l-l 

q 

E 

i=0 fe=0 
9+1 



q+2\ 



a 



(51) 



(52) 



— = 3/2 = vra d + na E E '*''*= cos(fcx)cr* 

i = l fe = l 

+ o(a^+^), 

where 

r 

y-- 

a 

and the parameter d introduced in (51) represents the post- 
Newtonian contribution to the radius of the ring to lead- 
ing order in a. The relevance of this parameter will be dis- 
cussed below. An important property of the Fourier series, 
is that they terminate due to the expansion in a. The func- 
tions Ii{y,x) form an axially symmetric set of solutions to 
Laplace's equation, regular everywhere except at r = and 
vanishing at infinity, see Petroff & Horatschek (2008). The 
set Ki{y,x) is the analogue for the generalized Laplace's 
equation (31), that can be found using separation of vari- 
ables in toroidal coordinates {rj,S^,tp). These are defined by 

b sinh ^ 



cosh ^ — cos rj ' 

6 sin 77 
cosh 5 — cos 77 ' 



(53) 
(54) 
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The set Ki can be defined by 

TT (cosh 5 -COS 77)^/2 1 
-^1 := ^^:t -^-r^ -P-i/2(coshg, (55) 



^263 sinh^ 



K, 



l_d 
bdb 



iCi 



with 



Ad d siny 9 

— +COSX: 



(56) 



(57) 



dfe 9& dr r dx 

On the axis, they have the simple form 

^'M^-^^lsTTT^ (on axis). (58) 

The idea now is to expand these functions with respect to 
a and then proceed to solve the field equations iteratively. 
We demonstrate the technique with the leading order in a. 



5 THE SOLUTION TO THE EQUATIONS 

The post-Newtonian correction to the surface to leading or- 
der in a is 



y2{x) = 7ra^[d-f o(a")]. 



(59) 



As with yo, this function is also independent of x to leading 
order and fulfils (38). Inside the ring, (32) becomes 



'^™«'+ 1^ + 16 = 0, 



(60) 
(61) 



dj/2 y dy 

whose regular solution is 

^00 = -~4:y'^ + ci 
The transition conditions are satisfied for 

ci = 8A + 4 (62) 

and where the constant Gio from (47) is 

Gio = -16. (63) 

The generalized Poisson equation (34) becomes 

4^ + i^ + 20y^-24A-58 = (64) 

dy^ y dy 

and we can immediately write down the regular solution 



Noo^-jy +(6A + yly +C2. 



Here, the transition conditions lead to 

C2 = -12A^ - 54A - 4Ad - 8d 

and for Fio from (45) 

Fio = -12A-4d-24. 



245 



(65) 



(66) 



(67) 



The equilibrium condition (36) gives an expression for 
the pressure coefficient goo that can be combined with (43) 
to give 

goo = ly* - (4A + 14)y' + 4A + ^ + 2d, (68) 

thus completing the solution to leading order in a. 

In general, stationary and axially symmetric fluids with 



a given equation of state are characterized by two param- 
eters, e.g. gravitational mass and angular momentum. The 
corresponding post-Newtonian expansions contain two ad- 
ditional free parameters, which arise because of the freedom 
one has in choosing which two quantities contain no post- 
Newtonian contribution, see Bardeen (1971) and discussion 
therein. We too would obtain two free parameters, but have 
already chosen to have b remain unchanged to first order in 
e and thus only have d as a free parameter in the solution.^ 
The leading order in a is equivalent to an infinite cylin- 
der, as in the Newtonian case, see Ostriker (1964a,b). In 
Bicak et al. (2004), relativistic, homogeneous cylinders and 
their post-Newtonian expansion are investigated. In that pa- 
per, there is no post-Newtonian contribution to the central 
pressure, which corresponds to choosing d = — 2A — 49/8 in 
our case. In order to compare our results to theirs, we intro- 
duce the proper radial coordinate and circumferential radial 
coordinate 



Vp 



Vc 



r y 

- / i/^dr' = / e" Ay' 


27r 
-^ / / J / ct 



(69) 
(70) 



(on surface) (71) 



The values of these coordinates on the cylinder's surface 
(using the choice for d given above) are found to be 

yp = l-|^^V[l + o(a")]+0(e*), 

yc = l-y^aV[l + o(a»)]+0(e*), 

in agreement with (6.10) and (6.11) in Bicak et al. (2004). 
The expressions for the pressure, e.g. as a function of proper 
radial coordinate, can also be shown to agree. 

The expansion to higher orders in a presents no particu- 
lar difficulties and we have carried it out up to seventh order. 
The resulting coefficients can be found in Appendix A up to 
fourth order, but the results presented in the next section 
are given up to seventh order. 



6 RESULTS 

In the following it will be convenient to use dimensionless 
quantities, denoted by a bar: 



Q^ 


1 


M 


Mb 


Q? 


G^' 


M ^ 


Mb 


J 

J ' 


gV 

C5 ' 


b 

b ^ 


g z 



Ql/2^1/2 



(72) 



(73) 



where M is the gravitational mass, Mb the baryonic mass 
and J the angular momentum. In particular 



6 = 



/Gjlb 



(74) 



holds. In Table 1, the Newtonian and first post-Newtonian 
results for the ring with the prescribed values gi/go = 0.7 



^ Note that to leading order in a, one only has the freedom to 
fix a single parameter. The parameter b does not appear to this 
order and we retain the freedom to choose d. 
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Table 1. Comparison between the Newtonian and first post-Newtonian (pN) results up to seventh order in a. The relative errors are 
computed by comparing with a highly accurate numerical solution with the prescribed values Qi/do = 0.7 and Zq = 0.05. For the 
definition of do,, dj and dM, see (76). 



Numerical value 






Relative error 






Relativistic 


Newtonian 


pN:d = 


pN: d = da 


pN: d = dj 


pN: d = dM 


n=4.9108 X 10~l 


-2 X 10-2 


-1 x 10-3 


-3 X 10-'' 


3 X 10-3 


2 X 10-3 


J=2.3168 X 10""' 


1 X IQ-i 


1 X 10-2 


8 X 10-3 


-8 X 10-3 


-2 X 10-3 


M=7.9661 X 10-3 


9x 10-2 


1 X 10-2 


7 X 10-3 


-5 X 10-3 


-1 X 10-3 


Mb=8.0842 X 10-3 


7x 10-2 


1 X 10-2 


5 X 10-3 


-4 X 10-3 


-1 X 10-3 



and Zq = 0.05 are compared with a highly accurate numeri- 
cal solution in full General Relativity. In Newtonian theory, 
there is (strictly speaking) no redshift, but the leading term, 
Zq — —V{)/(?, contains only the Newtonian quantity Vo and 
the speed of light, showing the relativistic origin. This equa- 
tion explains what is meant by the 'redshift' of a Newtonian 
configuration. The free parameter d is chosen in a number of 
different ways in the table. For one, we simply choose d = Q. 
Alternatively, we expand a chosen quantity 



J = Jn + JpN£ + 0{e ), 



(75) 



and define d so that there is no post-Newtonian contribu- 



tion 



ripN(d = dn) = 


= 0, 


JpN(d = dj) = 


= 0, 


V/"pN(d = dM) = 


= 0, 



(76) 



The solutions to these equations when expressed as series in 
a are 



dti = 



144A2 + 240A + 119 

24(4A + 1) 

6912A3 + 864A2 - 9684A - 5903 2 , 3^ 
+ 1152(4A + 1)2 " +''("') 

624A2 + 3120A + 2009 
" 24(12A + 7) 

202752A3 + 1487136A2 + 1765236A + 549479 2 

H rrr^ cr 



1152(12A + 7)2 



+ o(a3) 



1 ^ 

dM = -3(A + 2) + —a^ + o{a^) 
d4 



(77) 



Throughout Table 1, the post-Newtonian expressions 
yield a clear improvement compared to the Newtonian ones. 
The deviations from the accurate numerical values have two 
origins: Terms in a and terms in e are missing. Table 2 
demonstrates that for rings with Qi/Qo -^ 0.7, q = 7 is suffi- 
cient for relative errors smaller than lO"* in the Newtonian 
case. One can estimate that the error due to terminating the 
series in a at the order q results in a relative error of about 



((j| In a])', where it is convenient to use 



2 V eo 



(78) 



If the relativistic efi'ects are large in comparison (as is the 
case in Table 2 although Zq is only 0.05), then truncating 
in a is negligible. A rough estimate for e can be found by 
taking the leading term for the redshift in both e and a, 
inserting (78) and approximating slightly to find 



Zq 



1 — Oi/Qo V 6 - 2 ln(l — gi/go) 



(79) 



The relative factor between the Newtonian and post- 
Newtonian contribution to fl is about 5| In a\{a£)'^ ~ IO-2 
for the values from Table 2. 

Tables 3 and 4 present analogous results for Qi/Qo ~ 
0.95 and Zo — 10"3. The values of each expansion param- 
eter (i.e. a and e) are small enough that extremely high 
accuracy is observed. This also permits us to see how the 
results improve as one proceeds to the post-Newtonian level 
and with increasing orders in a. Tables 1 and 3 suggest that 
the choice d — dx is particularly accurate if the quantity x 
is prescribed. 

Figure 2 provides an example of the meridional cross- 
section of the surface shape of a ring. Both post-Newtonian 
results given there can be seen to constitute a marked im- 
provement as compared to the Newtonian one. Moreover, 
we find that the choice d = dM gives a significantly better 
result than d = 0, which is also the case for most values in 
Tables 1 and 3. 

A similar behaviour can be found in Fig. 3, where the 
pressure in the equatorial plane is plotted. Again the post- 
Newtonian results demonstrate an improvement as com- 
pared to the Newtonian, and here too the choice d — dM 
leads to better results than d = 0. 

The quality of the approximations can be seen for se- 
quences of rings with the prescribed radius ratios Qi/Qo = 
0.95 (Fig. 4) and Qi/Qo = 0.7 (Fig. 5). The relativistic curve 
'begins' in the Newtonian limit at the origin and 'ends' in an 
extreme Kerr Black Hole at the point (1, 1). The parametric 
transition of rings to Black Holes was first studied in Ansorg 
et al. (2003b) and had already been considered analytically 
for discs (Neugebauer & Meinel 1993; Meinel 2002). It is not 
surprising that the accuracy of the Newtonian and post- 
Newtonian curves in the vicinity of this ultra-relativistic 
point (at which Zq — > 00) is so poor. 3 In contrast, the highly 



2 Nonetheless even for this quantity the values in the Newtonian 
and post-Newtonian column are different, because Qi/Qo = 0.7 
and Zq = 0.05 are prescribed, and not the chosen quantity. 



3 It is a mere coincidence that some of the curves almost pass 
through the point (1, 1). 
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Table 2. Comparison between the Newtonian and first post-Newtonian (pN) results for the dimensionless angular velocity fl for different 
orders q in a. The prescribed parameters for the ring are Si/Qo = 0.7 and Zq = 0.05. The numerical calculations show that the value for 
Q for the relativistic ring is finum = 0.49108 (cf. Table 1) and for the Newtonian one is fJN.num = 0.48109 (see Table 1 in Ansorg et al. 
(2003b)). For the definition of dn, dj and djvfi see (76). 



g 


Newtonian 


pN;d = 


pN; d = dn 


pN: d = dj 


pN: d = dM 


1 


0.50085 


0.51126 


0.51192 


0.51364 


0.51301 


3 


0.48274 


0.49206 


0.49264 


0.49423 


0.49369 


5 


0.48135 


0.49059 


0.49116 


0.49277 


0.49223 


7 


0.48114 


0.49035 


0.49092 


0.49255 


0.49200 


20 


0.48109 


— 


— 


— 


— 



Table 3. Comparison between the Newtonian and first post-Newtonian (pN) results up to seventh order in a. The relative errors are 
computed by comparing with a highly accurate numerical solution with the prescribed values Qi/qo = 0.95 and Zq = 10^"'. For the 
definition of i^q, dj and dj^f , see (76). 



Numerical value 



Relative error 



Relativistic 


Newtonian 


pN:d = 


pN;d = dn 


pN: d = dj 


pN: d = dM 


n=9. 6226144 X 10"^ 


-5 X 10-4 


-1 X lO-'S 


-3 X 10-^ 


8 X 10-^ 


7 X 10-^ 


J=2.9191516xl0-'^ 


3 X 10--^ 


7 X lO-'' 


4x 10-« 


-2 X 10-*' 


-2 X 10-« 


M=8.5891100 X lO--"' 


1 X 10-^ 


5 X 10-" 


2 x 10-« 


-1 X 10-*' 


-7x 10-^ 


Mb=8.5914515 x lO-'""' 


1 X 10-^ 


4 X lO-** 


2 x 10-"^ 


-9 X 10-^ 


-8 X 10-^ 



0.08 



0.06 - 



Gfi^gl 



2 0.04 - 




Figure 3. The pressure profile of the ring with A = 0.7 and 
Zq =0.1 in the equatorial plane. For the definition of d^, dj and 
dM, see (76). 



1.2 
1.0 
0.8 
2riM 0.6 
0.4 
0.2 - 



Newtonian 




post-Newtonian {d - 


= 0) 


post-Newtonian {d - 


--dM) 


relativistic 






0.2 0.4 0.6 0.; 



1.0 1.2 



Figure 4. The sequence of rings with Qi/qo = 0.95 for g = 7. For 
the definition of dn, dj and dM , see (76). 



relativistic regime can be probed quite accurately using very 
high orders of the post-Newtonian approximation, at least 
for the disc of dust (Bardeen & Wagoner 1969, 1971; Petroff 
& Meinel 2001). In the Newtonian regime, M^/J -)■ 0, the 
curves become indistinguishable, as they must. Naturally, we 
require that there exist some neighbourhood of the Newto- 
nian limit in which the post-Newtonian approximation pro- 
vides an improvement compared with the Newtonian results. 
In order to check this, we have plotted the relative errors 
Ahm = 1 — (nM)/(r2M)reiativistic iu Figs 6 and 7. The errors 
corresponding to the sequence Qi/go = 0.95 shown in Fig. 6 
demonstrate the required behaviour. It is particularly ap- 
parent that d = dM produces better results than d = 0. For 
£'i/£'o ~ 0.7, the Newtonian limit shown in Fig. 7 seems to 
be in contradiction to the requirement: the post-Newtonian 
approximations in this plot provide no improvement. For 



this relatively small radius ratio, the truncation of the series 
at the order q = 7 provides the limiting value to the accu- 
racy of ~ 10-* as discussed above. On the other hand, the 
limitation in the relative error of ~ 10-* due to g = 7 for 
Qi/go = 0.95 is beyond the accuracy shown in Fig. 6. 

We have seen, that as with spheroidal bodies, the post- 
Newtonian approximation for toroidal ones can be used to 
great advantage. It is noteworthy that this work is possible 
without the existence of a known exact ring solution, even 
in Newtonian theory. It was therefore necessary to expand 
about the thin ring limit, which we were able to do up to high 
orders, making the results applicable even to fairly 'thick' 
rings. 
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Table 4. Comparison between the Newtonian and first post-Newtonian (pN) results for the dimensionless angular velocity Q for different 
orders q in a. The prescribed parameters for the ring are g\/go = 0.95 and Zq = 10^"^. The numerical calculations show that the value 
for n for the relativistic ring is finum = 0.096226144 (cf. Table 3) and for the Newtonian one is fSt-jnum = 0.096177081. For the definition 
of dn, dj and djv/i sec (76). 



Newtonian pN: d = Q pN: d = dn pN: d - 



pN; d = djyi 



1 0.096333800 0.096383025 0.096383103 0.096383208 0.096383200 

3 0.096177653 0.096226613 0.096226690 0.096226795 0.096226786 

5 0.096177085 0.096226044 0.096226121 0.096226226 0.096226217 

7 0.096177081 0.096226040 0.096226117 0.096226222 0.096226214 



20 0.096177081 





Figure 2. The cross-section of the ring with M = 0.05 and Zq = 0.2. The solid line is the highly accurate numerical result. The ' 
dots, which coincide with it almost exactly, represent the post-Newtonian result with d = dM (cf. (76)). The thin dotted line is the 
post-Newtonian result with d = and the dashed line is the Newtonian result. All analytic curves were generated with q = 7. 




1.0 0.2 0.4 0.6 0.. 



1.0 1.2 



Figure 5. The sequence of rings with Qi/ Qo = 0.7 for g = 7. The 
line types are the same as in Fig. 4. 
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APPENDIX A: THE COEFFICIENTS 

In the Tables A1-A6 the post-Newtonian coefficients up to 
the order g = 4 are presented. 
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Table Al. Coefficients Lj up to the order g = 4. 

i Li 

1 __;i== (i44A2 + 96Ad + 240A + 24d+119) 

2 

3 --rr7TT77TTTT372 (3744Ad + 1728A2d + 1439 + 7236A + 6048A^ + 2076d + 11448A2) 

4 

Table A2. Coefficients K^fe up to the order q = A. 

t K l^ik 

1 1 

2 1 

2 2 |A^ + fdA+^A+^ + id 

3 1 

3 2 

O O 32 "■" 32 3072 73728 768 

4 1 

, r. 25\3,39\2,25j\2 515 \ , 315 ,j\ 2093 1115 J 

* ^ "8"-^ + 16-^ + — C(A - 1152'^+ 128"-^^ 6912 ^ 92T6" 

4 3 

4 4 75 \3 I 8111 ^2 , 375 J^2 , 2725 ^\ , 328073 ^ i 4190617 , 5885 j 
32 ^ "■" 2304 "■" 256 ~'~ 2304 ^ ^ 184320 ' 13271040 ' 55296 
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Table A3. Coefficients N^i_{y) up to the order q = A. 



3 



4 



iV' 



ik 



-12A2-54A-4dA-245 _8rf+(6A + f)y2_ |y4 

1 

1 1 (-6A2-16A-2dA-f -d)2y+(iA+|)j/3-^j/5 



2 |A^ + idA+liA + |d+i^ + (-iA^-6A-idA-id-^)y^+(i|A+^)jy^-^I/' 

2 1 

2 2 (-^A^-2dA-fA-iJ^-^d)y^+(iA+^),^-^,« 



85 -v 2 I 13 j^ I 343 -v , 11987 , 91 A,, , ( 13 a, 2 7>a 1727 x 25 j 24167 -) ,,3 , 

221 \ -L 1307 '! ,,5 25_ 7 

192 ' 576 / y 128 y 



r-22\2 _ 15 ji, _ 14185 X _ 50857 , _5_j-, „3 , /' 515 ), , 16465 -) „5 _ X,,^ 
V 96 32 "^ 4608 ^ 36864 "■" 768 "/ » "•" V 768 ^ 18432 / « 96 y 

125 \4 _ 75^1,3 _ 1745 1,3 _ 47291 1,2 _ 55ji,2 _ 1097 j-i _ 623875 i, , 1633 ^ _ 118327 , 
16 '^ 32 "'^ 64 ^ 1536 ^ 16 ^ 1536"'^ 55296 ^ ' 9216 221184 ' 

(11 \3 _l_ 405 1,2 I 3767 i, _, 13 J\ _i_ .9]_ J _, 9767 'i 2 , (_39 \2 _ 873 \ _ 21_ j\ _ _25_ j _ 26059 'i 4 _,_ 

164^^ """256^ "•" 3072 '^ "•" 32 "^ "•" 384 ""•" 12288 / » "*" I 32 '^ 256^ 64"^ 256 12288 / » "*" 

(-135;^, 1415 -) 6 _ .675. 8 

>,256 ^ 1536/*' 8192 W 



_ 55 \3 _ 5 j\2 _ 3517 1,2 _ 151 ja _ 246605 i, , 534997 , 5329 ^\ ,2 , 
gA 2"^ 256^ 128 "'^ 36864 ^ "•" 884736 "•" 9216 "/ y "'" 

— 211x2 _ 185 J\ _ 6933 1, _ 985 j _ 1142287 \ 4 , 1- 2389 1, , 95237 -1 6 _ 175 8 
384 '^ 384 "-^ 2048 ^ 9216 " 442368 J ^ ' \ 3072 ^ "•" 73728 / » 1536 » 

) 

_ 2821 1,2 _ 455 j\ _ 418169 \ , 1001 j _ 10764521 \ 4 , 1- 2611 \ , 85229 ^i 6 _ 133 8 
2304 2304 184320 ~'~ 55296 13271040/" """ 1 9216 ~'~ 221184 / » 4096 » 
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Table A4. Coefficients Fn up to the order q = A. 

I i Fii 

1 -12A-24-4d 

1 1 

1 2 ^ 
^ ^ 16 

1 3 

-I , 125 A 3 75i\2 555 ^ 2 17357 \ 125 7\ 115483 175 j 
^ 4 16"'^ ^32"-^ ^32"^ ^ ^536"-^ ^ "eT"-^ ^ 55296 ^ 5T2" 

2 1 fA+f+d 

2 2 

135 \2 I 1939 



o q 135 -V 2 , 1939 \ I 76^\ , 125^ , 1M527 

2 -i -g-A + -(j^A+ — fflA+ ^J2-"+ T536" 

2 4 

o r, 25 \ 2 15 7-, 23-1 3 J 2693 

i 2 - — A - — dA- — A- jCf- ^76" 

3 3 

o . __ 55 \3 _ 4375 \2 _ c,J\2 _ 168149 \ _ W21t]\ „ 3024601 _ 2645 j 

■^ ^ 8 256 "-^ 12288 192 "^ 884736 2304" 

A ■^ 143x2 . 5 ,., I 11729 N , 4801 , 23. 

t o 192 '^ T" 16 "'^ ^ 9216 '^ ^ 8192 ^ 96" 

4 4 

c- 4 —lis \3 _ 13427 a2 _ 125 ja,2 _ 2125 j\ _ 1505639 -v _ 24750911 _ 44165 1 
192 6912 192 3456 1105920 79626240 331776 



Table A5. Coefficients wj;.(j/) up to the order q = A. 



i k wl). 



-4j/^ + 4 + 8A 

1 

1 1 (12A-4)jy-3jy3 

2 |A+f + (9A-3)y2-fj/^ 
2 1 

2 2 (lOA - i) y^ - |y^ 

3 

3 1 (-|A+f),+ (fA-i)y«-^,^ 

3 2 

o o fa^x 3395 \ 3 63 5 

■^ •^ I 48 '^~ 1152^2^ ~ 64?^ 

. n 25 n3 , 115 n2 I 1057 \ 1 115 , / 9 1, , 165 ^ „2 , ^375 x 475 \ 4 245 6 

4 ygA + -gg-A + 2304^+ TT52 + 1^32^+ "eTJ J' + 1^2"-^^ 128 j!' ^ 128^ 

4 1 

4 9 1^5x2 I 365 X , 30815 \ 2 , ^ 945 x 2905\ 4 315 6 

4 3 

4 4 flMS A - 22223\ 4 _ _693. 6 
V 576 13824/ » 1280'' 
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Table A6. Coefficients Gn up to the order q = A. 



Gh 








-16 






1 









2 


-12 






3 









4 


25 A 2 895 X 


11305 
1152 


2 


1 


6 




2 


2 







2 


3 


55 A I 409 




2 


4 







3 


2 


-5A-f 




3 


3 







3 


4 


rx2 995 X 


34361 
4608 


4 


3 


55 1, 1 869 
48^ ^ 1152 




4 


4 







5 


4 


25 1,2 2195 \ 
48 ■^ 3456 ^ 


20819 
103680 
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